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^ . Abstract 

Two-particle angular correlation for charged particles emitted in Au+Au 
collisions at the center-of-mass of 200 MeV measured at RHIC energies re- 
vealed novel structures commonly referred to as a near-side ridge. The ridge 
phenomenon in relativistic A+A collisions is rooted probably in the initial 
conditions of the thermal evolution of the system. In this study we ana- 
lyze the evolution of the bumping transverse structure of the energy density 
q I distribution caused by fluctuations of the initial density distributions that 

could lead to the ridge structures. We suppose that at very initial stage of 
collisions the typical one-event structure of the initial energy density profile 
can be presented as the set of longitudinal tubes, which are boostinvariant 
>• ! in some space-rapidity region and are rather thin. These tubes have very 

high energy density comparing to smooth background density distribution. 
The transverse velocity and energy density profiles at different times of the 
evolution till the chemical freeze-out (at the temperature T=165 MeV) will 
be reached by the system are calculated for sundry initial scenarios. 
PACS number(s): 25.75.Gz, 24.10.Nz 



X ■ 1 Introduction 

.S3. 

Measurements at the Relativistic Heavy Ion Collider (RHIC, Brookhaven) 
have revealed that the long range structure of two-particle angular correlation 
functions is significantly modified by the presence of the hot and dense matter 
formed in relativistic heavy ion collisions [I]. Novel correlation structures 
over large pseudorapidity interval Ay were observed in azimuthal correlations 
for the intermediate particle transverse momenta pr ~ 1 — 5 GeV/c [2J |3]. 
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First the striking "ridge" events were revealed in studies of the near-side 
spectrum of correlated pairs of hadrons by the STAR collaboration [I] . The 
spectrum of correlated pairs on the near side (defined by the trigger particle 
direction) extends across the entire detector acceptance in pseudorapidity 
interval of order At] ~ 2 units and is strongly collimated for azimuthal angles. 
Here At] is the difference in pseudo-rapidity 77 = — In (tan (6/2)), where 9 is 
the polar angle relative to the beam axis, between two particles and Atp 
is the difference in their azimuthal angle (p (in radians). The analyzes of 
measurements by the PHENIX [5] and PHOBOS [6] collaborations confirmed 
the STAR results. 

The discovery of the ridge has aggravated quantitative theoretical ana- 
lyzes which propose rather different explanations [UIH1IH]- The first, [7] treats 
the ridge as an initial-state effect. The authors point out that the enhanced 
two-particle correlations are a natural consequence of the correlations in the 
classical color fields responsible for multi-particle production in relativistic 
heavy-ion collisions. Due to fluctuations of color charges in colliding nuclei 
the longitudinally boost-invariant and transversally inhomogeneous structure 
of the matter can be formed. When it expends hydrodynamically it could 
lead to the ridges [8]. The others explore a final-state effect as the origin 
of the ridge [9]. The evolution of the system and other later-stage effects can 
modify these correlations, which, in fact, are associated with fluctuations in 
the energy densities at the final stage. The recent measurements by the CMS 
Collaboration [TU] observed unexpected effect in proton-proton collisions at 
the LHC. The clear and significant "ridge" -structure emerges at Aip m ex- 
tending to \Ar]\ of at least 4 units. This novel feature of the data has never 
been seen in two-particle correlation functions in pp or pp collisions before. 
This unusually elongate structure the ridge -remains after removal of elliptic 
flow and ordinary jet correlations [UJ. This is faced with the problem of 
causality which, probably can be solved only if one supposes that the ridge 
phenomenon in relativistic A+A collisions is rooted in the initial conditions 
of the thermal evolution of the system. The aim of this study is to check 
this hypothesis by an analysis of the developing energy density in the system 
which at very initial stage of collisions has transversally bumping tube-like 
fluctuations with boost-invariant homogeneous structure within some space- 
rapidity region. 
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2 Calculation of energy density distributions 



The numerical results presented in this section were obtained on the basis 
of original 3D ideal hydro-code, described in details in [12]. The analysis is 
based on hydrodynamic approach to A+A collisions and considered within 
the Boltzmann equations. It is consistent with conservations laws and ac- 
counts for the opacity effects. The hydrodynamic evolution starts at the 
time To- We use Bjorken-type initial conditions at t$: boost-invariance of 
the system in longitudinal direction, initial longitudinal flow vl = z/t with- 
out transverse collective expansion. In present calculation we compare the 
transverse-velocity profile of hydrodynamic flow and energy density profile 
which evolve in time till the chemical freeze-out (T = 165 MeV) in different 
initial scenarios. The one of them corresponds to the smooth Gaussian pro- 
file with radius R and energy density as it was considered in [13] at r = 0.2 
fm/c. The other scenarios are based on transversally bumping tube-like ini- 
tial conditions at r . These tubes are rather thin transversally and relatively 
long in the direction of beam axis; with radii a; = 1 fm. The general energy 
density distribution at tq could be written as: 

2 + 2 N t __(x-xi) 2 +(y-Vi) 2 

E = E b e~^ +Y,Eie =? , (1) 

i=Q 

Ri = a* = vl (2) 

where E b is the maximum of average energy density distribution, Ei are 
the maxima of tube-like fluctuations, Ri are the positions of the fluctuation 
locations and N t is the number of tubes. Instead of a study of the result 
over very many fluctuations, that should be finally averaged over azimuthal 
angular (it brings symmetry), we will be based here on the possible typical, 
or representative, initial fluctuation which are already maximally symmetric 
in azimuthal plane. The following initial configurations are considered: 

i The configuration without fluctuation: distribution of initial energy 
density corresponds to the Gauss distribution with R = 5.4 fm and 
maximum energy density at r = is Eb = 90 GeV/fm 3 ; 

ii The configuration with one tube (fluctuation) in the center: energy 
density profile is the Gauss distribution with a = 1.0 fm and the maxima 
value 270 GeV/fm 3 ; 
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iii The configuration with one tube (fluctuation) shifted from the center: 
E b = 90 Gev/fm 3 ; R = 5.4 fm; E = 270 GeV/fm 3 ; R = 3 fm; a = 1 
fm; The results of calculations are presented for initial time and for r 
= 1, 2, 10 fm/c ; 

iv The configuration with four tubes (fluctuations): Eb = 85 GeV/fm 3 ; R 
= 5.4 fm; Ei= 250 GeV/fm 3 ; Ri = 5.6 fm; Oj = 1 fm; The evolution of 
energy density profiles is presented on the Fig. 1; 




Figure 1: 3D plots of energy density profiles with 4 tube-like initial conditions 
for r = 0.2, 1, 2, 3, 6 and 10 fm/c 
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v The configuration with ten tubes (fluctuations): Ef, = 25 GeV/fm 3 ; R 
= 5.4 fm; #0 = fm; R { (i < 3) = 2,8 fm; R { (i > 3) = 4.7 fm; a { = 1 
fm; Et = 4E b exp(-Rl/R 2 ) . 

For all the considered cases the traces of the initial fluctuations - bump- 
ing final energy distributions - remain after the system evolution that should 
lead to a non-trivial structure in observed correlations. Besides the energy 
density profiles the evolution of transverse velocity profiles of hydrodynamic 
flow in different scenarios also was considered. The fluctuations of the initial 
conditions result in the fluctuation of the transverse velocity in the system. 
At the supposed thermalization time r = 1 fm/c the corresponding fluctua- 
tions in the transverse velocity averaged over azimuthal angular and radius 
are approximately 30 % between the cases when the high density fluctuation 
peak is just in the center and when it is shifted while at the later times 
(r = 10 fm/c) it is only 2-3 % . This may lead to the important conclu- 
sions as for small fluctuations of the mean transverse momenta of observed 
particles, despite big fluctuations of transverse velocities at thermalization 
time. 

3 Conclusions 

The basic hydrokinetic code, proposed in [T2] was modified to include the 
tube-like initial conditions with the aim to study how the initial correlations 
in the energy density evolve with time. We found that the effect of the 
fluctuations of the initial conditions does not wash out during the system 
expansion that, probably, leads to the ridges structures of the correlations, 
which are caused by these fluctuations. The evolution of the energy density 
distributions and velocity profiles are calculated in the framework of the 
3D hydrodynamics. Possible physically grounded configurations of initial 
density profiles are proposed. The first spectra calculations were done in the 
frameworks of the hydrokinetic model (HKM) [121 E! which allows describing 
all the stages of the system evolution as well as a formation of the particle 
momentum at the decoupling stage. The further investigations on this matter 
and description of correlations could be done in frames of this approach. 

4 Acknowledgments 

Yu.S. gives thanks to P. Braun-Munzinger for support this study within 
EMMI/GSI organizations. The researches were carried out in part within 



5 



the scope of the EUREA: European Ultra Relativistic Energies Agreement 
(European Research Group GDRE: Heavy ions at ultrarelativistic energies) 
and is supported by the State Fund for Fundamental Researches of Ukraine 
(Agreements No. F33/461-2009, 2011) and National Academy of Sciences of 
Ukraine (Agreements No F06-2010, 2011). The numerical calculations were 
done by using the cluster and GRID environment of the Bogolyubov Institute 
for Theoretical Physics of National Academy of Sciences of Ukraine. 

References 

[1] Alver B. et al. (The PHOBOS Collaboration) Phys. Rev. C 81 024904 
(2010). 

[2] Alver B. et al. (The PHOBOS Collaboration) Phys. Rev. Lett. 105 
022301 (2010). 

[3] Abelev B.I. et al. (The STAR Collaboration) Phys. Rev. Lett. 104 
062301 (2010). 

[4] Horner M. G. (for the STAR Collaboration) J. Phys. G: Nucl. Part. Phys. 
34 S995 (2007). 

[5] Adare A. et al. (PHENIX Collaboration) Phys. Rev. C 78 014901 (2008). 

[6] Wosiek B. (for the PHOBOS Collaboration) J. Phys.G 35 104005. 

(2008) . 

[7] Dumitru A. et al. Nucl. Phys. A 810 91 (2008). 

[8] Hama Y. et al. larXiv:0911.0811l [hep-ph](2009). 

[9] Schenke B., et al. J. Phys. G: Nucl. Part. Phys. 35 104109 (2008). 

[10] Khachatryan V. et al. (The CMS Collaboration) CMS-QCD-10-002, 
CERN-PH-EP/2010-031(2010). 

[11] Abelev B.I. et al. (The STAR Collaboration) Phys. Rev. C 80 064912 

(2009) . 

[12] Akkelin S.V., Hama Y., Karpenko Iu.A., Sinyukov Yu.M. Phys. Rev. 
C 78 034906 (2008). 



6 



[13] Sinyukov Yu.M., Karpenko Iu.A. and Nazarenko. A. V. J. Phys. G: Nucl. 
Part. Phys. 35 104071 (2008). 

[14] Sinyukov Yu.M., Akkelin S. V. and Hama Y. Phys. Rev. Lett. 89 052301 
(2002). 



7 



